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ABSTRACT 

We show that, by studying the arrival times of radio pulses from highly-magnetized pulsars, it may be possi- 
ble to detect light spin-0 bosons (such as axions and axion-like particles) with a much greater sensitivity, over 
a broad particle mass range than is currently reachable by terrestrial experiments and indirect astrophysical 
bounds. In particular, we study the effect of splitting of photon-boson beams under intense magnetic field gra- 
dients in magnetars and show that radio pulses (at meter wavelengths) may be split and shift by a discernible 
phase down to a photon-boson coupling constant of g ~ 10 _14 GeV _1 ; i.e., about four orders of magnitude 
lower than current upper limits on g. The effect increases linearly with photon wavelength with split pulses 
having equal fluxes and similar polarizations. These properties make the identification of beam-splitting and 
beam deflection effects straightforward with currently available data. Better understanding of radio emission 
from magnetars is, however, required to confidently exclude regions in the parameter space when such effects 
are not observed. 

Subject headings: elementary particles - pulsars: general - magnetic fields 
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1. INTRODUCTION 

The Peccei-Quinn mechanism (Peccei & Quinn 1977) was 
devised to elegantly solve to the strong-CP problem of QCD. 
This was accomplished by postulating a new quantum field 
and a new class of particles associated with it. The particles 
are spin-0 pseudo-scalars that couple very weakly to the elec- 
tromagnetic (EM) field. It later became apparent that such 
particles could also provide a solution to the dark matter prob- 
lem. To date, however, there is no observational evidence for 
the existence of such particles, dubbed axions, and it is not 
clear that the Peccei-Quinn solution actually works. 

Besides QCD axions there are also the putative axion-like 
particles (ALPs). Such class of hypothetical particles has re- 
lations to dark energy and may be related to the quintessence 
field. These particles have been proposed as a possible so- 
lution to the apparently low opacity of the Universe to hard 
y-ray radiation (e.g., de-Angelis et al. 2008) alongside other 
explanations (Reimer 2007, Jacob & Piran 2008). 

There is a longstanding interest in determining the physi- 
cal properties of axions/ ALPs. At present, laboratory experi- 
ments and astrophysical bounds imply that their coupling con- 
stant g < 1 0~ 10 GeV ~ 1 . Mass limits are less stringent: if QCD 
axions are concerned, then their mass is probably > 10~ 6 eV 
since otherwise the Universe would over-close, in contrast to 
observations. These limits, however, do not apply for ALPs. 

In recent years it has been realized that the Universe, pro- 
viding us with extreme conditions, can help us constrain the 
axion and ALP properties. For example, the implied trans- 
parency of the universe to y-ray radiation, if due to ALPs, 
puts constraints on the coupling constant and the mass. More 
recently, Chelouche et al. (2009; hereafter C09), have sug- 
gested a new way to detect axions by looking at the spectra of 
compact astrophysical objects such as magnetars, pulsars, and 
quasars. Unlike terrestrial experiments where the magnetic 
fields are relatively small and the system size is limited by 
human capabilities, magnetars have magnetic fields, B, which 
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are ~10 orders of magnitude stronger than laboratory ones 
and extend over R ~ 10km scales rather than a few meters. 
As the photon-particle conversion probability is oc B 2 R 2 , more 
stringent constraints on axion/ ALP physics may be obtained. 

In this paper we propose a new method to detect ALPs 
by utilizing the photon- ALP duality (Guendelman 2008a,b,c). 
As we shall show, if ALPs exist then it is possible to de- 
tect their signature in the light-curve of magnetars down to 
g ~ 10 _14 GeV _1 ; i.e., much more sensitive than other means 
and complementary to the spectroscopic method of C09. This 
paper is organized as follows: in §2 we lay out the formal- 
ism used to calculate first order effects in g resulting in beam 
splitting effect that are analogous to the Stern-Gerlach exper- 
iment. Order of magnitude estimates are given demonstrating 
the feasibility of the method for detecting ALPs. Section 3 
considers a more quantitative model pertaining to magnetars, 
and shows the predicted light curves when photon- ALP cou- 
pling is important. We conclude in §4. 

2. SPLITTING IN IN-HOMOGENOUS MAGNETIC FIELDS 

In a recent series of papers, Guendelman (2008a,c) pointed 
out a photon-particle duality in a system containing a photon- 
particle interaction term of the form 



^gF^F^a^gE-Ba 
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where E is the electric field (associated with the photon), B 
the magnetic field, and a the axion field, g is the coupling of 
particles to the EM field. The full Lagrangian for the system 
can be written as 
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L = —F^F^ - -m 2 ^ 2 + -d^ad^a - -m 2 a a 2 +L iat (2) 

which is the free EM Lagrangian including an effective mass 
term which takes into account potential refractive index in the 
medium, as well as the Klein-Gordon equations for free par- 
ticles (m a ,my are the particle and effective photon mass, re- 
spectively). In the absence of L mt , photons and particles (e.g., 
axions) are well defined energy states of the system. However, 
once the interaction term is introduced, the only well defined 
energy states of the system are mixed photon-particle states. 
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FIG. 1. — Photon-boson beams of different "charges" (see Appendix and 
Eq. [5) would be split along magnetic field gradients in a way similar to the 
Stem-Gerlach experiment. Photons propagate along the z-axis with their po- 
larization along magnetic field lines (x-axis). A Schematic view of a magne- 
tar is also shown. Magnetic field lines originate from the magnetic pole with 
plasma in its vicinity emitting beamed radiation. One possible orientation of 
the coordinate system is also shown. 

The equation of motion for the photon-particle system takes 
the form (e.g., C09, Raffelt & Stodolsky 1988), 3 



k 2 - co 2 
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where co is the photon energy and Bii the magnetic field in 
the direction of the photon polarization (the photon's E field). 
Clearly, neither pure photon nor pure ALP states are eigen- 
states of the system but rather some combination of them. 
Let us now focus on the limit 

\m 2 a -m y {co) 2 \ «£B||0)~ KTVwfiieA^eV 2 (4) 

where Bi 6 =B||/10 16 G,g= lO^V^GeV- 1 , and the pho- 
ton wavelength, X = A m m. This condition is met either near 
resonance where m 2 ~ m 2 or when both masses are individu- 



ally smaller than y/gB\\ co (which limit is actually met is irrel- 
evant). The eigenstates of equation [3] are then given by 
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where \a) is the axion state and |y) is the photon state. The 
eigenvalues are ra 2 = ±gB\\CB. By analogy with optics, these 
masses are related to effective refractive indices: n± = 1 + 
8n± ~ 1 — ot 2 /2(0 2 (for |5n±| -C 1) meaning that different 
paths through a refractive medium would be taken by the rays. 
We note that there is no dependence on the particle or photon 
mass so long as equation|4]is satisfied. 

In terms of the refractive index, the equation of motion for 
a ray may be found by minimizing the action / dsn(s). This is 
completely analogous to mechanics where we substitute £ — * 
COn±. In this case, a force is <9£/<9s = ±(g/2)(dB/ds). Using 
our simplified geometry, depicted in figure 1 , the momentum 
imparted on each state is 



Spt = T(g/2) dz(dB x /dy) 



(6) 



where B x = B x (y,z) (note that dt = dz in the adopted units). 
Clearly, each of the beams will be affected in a similar way 

3 We work in natural units so that % = c = I. 



while gaining opposite momenta so that the total momen- 
tum is zero and the classical wave packet travels in a straight 
line (along the z-axis). This effect is analogous to the Stern- 
Gerlach experiment (see Fig. 1). Beam splitting effects, aris- 
ing due to photon-particle (polariton) mixing, have been mea- 
sured in the laboratory (Karpa & Weitz 2006) and that an anal- 
ogy exists between this case and scalar QED (see appendix). 
In the limit n± ~ 1, the separation angle between the beams 



e~2p- i \8p y \~co- [ gf G B 
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where p is the beam momentum along the propagation direc- 
tion, i.e., the z-axis. This expression holds for small splitting 
angles and assumes relativistic axions. We also approximated 
/ dz(dB x /dy) = f G Bu where f G (N) is a geometrical factor 
depending on the magnetic field geometry, the inclination of 
our line-of-sight through the magnetized region (e.g., for pul- 
sars and magnetars the magnetic field is predominantly dipo- 
lar and fa < 1; see §3), and on the photon polarization. An 
implicit assumption in the above expression is that the field is 
monotonically increasing or decreasing. In situations where 
the magnetic field is stochastic, 8 p y can no longer be evalu- 
ated according to equation [6] which is linear with distance (or 
time) and a better treatment is that of a random walk nature 
whereby the average 9, (9) \/t; such cases are beyond the 
scope of this paper and are likely to be less relevant to com- 
pact astrophysical objects whose magnetic fields are thought 
to be relatively ordered. 

Assuming splitting angles 9 ~ 10~ 2 6_2rad are detectable 
at radio wavelengths (see §3), then the minimum coupling 
constant which can be probed, provided equation|4]holds, is 



g min ~2x 10- 14 A m V G Xe^GeV 



(8) 



It is important to establish that the above condition on the 
photon and axion mass can be materialized under realistic 
conditions. Consider the non-resonant case: as the mass of 
ALPs is unknown, there is no reason to suspect such a con- 
dition is irrelevant. As for the photon mass, there are two 
contributions to the refractive index of the photon in magne- 
tized plasma: vacuum birefringence and a plasma term (e.g., 
C09). The condition for the plasma term to be negligible is 



f G 9- l {co p /cof^\. 



(9) 



The fact that we observe a signal at photon energies CO, means 
that CO » CO,,. Taking f G ~ 0.1 and 9 ~ 10~ 2 rad, this condi- 
tion is met if Co/co p > 5 (see §3). We shall therefore neglect 
plasma effects in our analysis. A more restrictive condition 
is associated with the vacuum birefringence term which was 
calculated by Adler (1971) and used in the context of photon- 
particle oscillations by C09. In particular, if the refractive 
index n = 1 + 8n9 ED then we demand that 



20- 1 5« QED / G <l. 



(10) 



For 9 = 10~ 2 rad, f G = 0.1, and <5n QED ~ 0.05 (C09) we get 
that the above ratio is of of order unity. While this condition 
is not strictly satisfied, we note that for slightly larger 9 val- 
ues (e.g., slightly lower co), this condition is met. The proper 
treatment of cases in which the mass terms are of order the 
mixing terms is beyond the scope of this paper and is the sub- 
ject of a follow up investigation. In resonance, the QED and 
plasma contributions to the effective mass cancel; we elabo- 
rate more on these conditions in S3. 
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FIG. 2. — Splitting of the peak pulse at = rad in a single rotation light- 
curve of XTEJ1810197 (dotted line; see Camilo et al. 2006) for several 
values of the coupling constant, g and assuming A = 1 m, /g^II = l"' 5 G. 
Note the similar fluxes of the split signals whose sum corresponds to that of 
the original pulse. Looking for the effects of pulse-splitting in the radio light- 
curves of magnetars allows one to be considerably more sensitive to light 
bosons compared to e.g., the CAST experiment and other astrophysical con- 
straints. In particular, pulse splitting at meter wavelengths can be detected 
down to coupling constants g m \ n > 10" 14 GeV"' form n < 10~ 7 eV. Observ- 
ing at longer wavelengths (and assuming all other parameters are fixed) will 
proportionally increase the sensitivity to lower values of g mm (see Eq.[8j. 

For splitting to be observed, we require that the photon- 
boson beams be able to propagate through the medium hence 
that the refractive index is real, i.e., 

ml/(0 2 =gB\\/(0~e/f G < 1, (11) 

Clearly, for 9 > 0.1 rad (taking fc = 0.1), one beam would 
be attenuated as it travels through the medium. In such cases, 
it is more appropriate to talk about pulse shifting, as only 
one beam is deflected from its original path and remains ob- 
servable. In this case the deflection angle, or phase-shift, is 
9 = 9/2. 

The time-delay between the two split beams is St ~R\m+ — 
m 2 _\/2a> 2 c<s:R/c~ 10~ 5 (/?/10km) s. Such small splitting is 
unlikely to be observable given the typical duration of pulses 
from magnetars (see below). 

3. ALPS AND THE RADIO LIGHT-CURVES OF MAGNETARS 

Pulsed radio emission from magnetars, whose rotation pe- 
riod is of order seconds, was discovered only very recently 
(Camilo et al. 2006) showing numerous narrow pulses over 
individual rotations. The phase duration of individual pulses, 
8(j) ~ 10~ 2 rad and so they become visible in only a small 
fraction of the period, 8^ ~ 10~ 3 (see e.g., Camilo et al. 2006 
for the case of XTE J1810197). Given their high luminos- 
ity, the pulsed emission is likely to be beamed (with beaming 
factors of order > 100). The radio emission is known to be 
linearly polarized (Camilo et al. 2006). 

Wishing to gain better understanding of the phenomenology 
of beam splitting, we note that there is no complete model for 
the radio emission from magnetars, and that the magnetic field 
configuration is not well known. Due to the inherent uncer- 
tainties we shall consider a very simple picture of a magnetar 
and neglect the, possibly important, effects of magnetic loops 
(Thompson & Duncan 2001) and global field twists (Thomp- 
son et al. 2002). In our model, the radio emission is beamed 
and originates from regions with high magnetic fields. The 
magnetic field is assumed to be dipolar and the rotation axis 
misaligned with the magnetic axis (see Fig. [TJ. 

We have numerically integrated equation [6] from the emis- 
sion point of the photon near the surface of the neutron star 
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up to large radii assuming various inclination angles of the 
magnetosphere with respect to the line-of-sight (photon was 
assumed to be linearly polarized at some direction). Perform- 
ing various such integrations we choose to adopt fc ~ 0.1 
and show the effect of splitting on sub-pulses in figure [2] 
Clearly, the effect is large and could be easily observed. Each 
pulse is split into two sub-pulses with equal fluxes (since 

KtIvMI 2 = l(y|V-)| 2 when equation 0] is fulfilled) and hav- 
ing a discernible phase difference. At the CAST limit (g < 
10 _10 GeV _1 ; Andriamonje et al. 2007), phase shifts would 
be of order unity. Most notably, even for g > 10~ 14 GeV~', 
the effect can be measured at >meter wavelengths. Whether 
or not non-split pulse components would be visible, depends 
on the photon polarization: splitting does not affect photons 
whose polarization is at right angles to the magnetic field di- 
rection. As splitting increases, asymmetric splitting, with re- 
spect to the original non-split pulse, may arise. This is caused 
by variations in the magnetic field gradients along different 
sight-lines through the magnetosphere and, in case |m^|/a) 2 
is of order unity, also the somewhat different propagation 
speed of each beam through the medium. 

The effect of splitting is wavelength dependent. In particu- 
lar, 9 oe X and so, by comparing the light-curve of magnetars 
at multiple bands, one may be able to identify split pulses 
from double pulses (i.e., those which are emitted as such at 
the source). This naturally assumes that the emission region 
for all bands is similar. For large enough splitting angles (and 
depending on fo), one beam will be attenuated and only pulse 
shifting by a phase 9 can be measured. To quantify the phase 
shift one needs to either measure it with respect to the origi- 
nal pulse (e.g., if some photons with a different polarization 
were not deflected by the effect), or statistically, by measuring 
typical pulse shifts between various bands and identifying the 
proper wavelength dependence (assuming there are no intrin- 
sic systematic phase shifts of the pulses between the various 
bands). 

The formalism adopted here holds equally well for the case 
of resonances whose importance is in the fact that they al- 
low, in principle, to probe more massive ALPs/axions. As 
the photon mass is wavelength dependent, resonance would 
occur at particular frequencies, Xq, over a limited spectral 
band. Once resonance bands are determined (depending on 
the plasma density in the magnetosphere and magnetic field 
strength; C09), beam splitting effects can be evaluated us- 
ing equation [8] with X = Xq. As discussed in C09, typical 
wavelengths for resonances which would probe more mas- 
sive ALPs/axions are in the sub-mm to infrared range hence 
the expected splitting angles are several orders of magnitude 
smaller (assuming all other parameters are fixed) and the sen- 
sitivity to low values of g is correspondingly lower. 

3.1. Caveats 

While we have demonstrated that the presence of finite 
photon- ALP coupling could lead to an easily detectable signa- 
ture in the light curves of magnetars, failing to detect splitting 
features cannot be used to place reliable limits on g so long 
as the radio emission mechanism in magnetars is not well un- 
derstood. For example, it might be that magnetic fields in the 
vicinity of the radio emitting region are considerably smaller 
than assumed here (e.g., if photons are emitted high above 
the stellar surface, in the outer parts of the magnetosphere). 
It could also mean that the photon polarization is less favor- 
ably inclined with respect to the magnetic field. Lastly, it may 
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FIG. 3. — The axion/ALP parameter space which can be explored by beam 
splitting in the radio light curves of magnetars as compared with other de- 
tection methods for several values of £ = (/gB||/10' 6 G)(A/m)]. Clearly, 
beam splitting/shifting is a potentially very efficient and sensitive method for 
detecting ALPs at the low mass end where equation[4]is satisfied (resonances 
are not shown; see text). In particular, the sensitivities gained by this method 
are 3-4 orders of magnitude better than those reached by CAST and by dating 
of horizontal branch (HB) stars. Dashed thick lines mark the g-value above 
which one beam attenuation occurs (for each f and assuming fo =0.1). 

be that the radio emission is less beamed than implied by the 
pulse duration, in which case splitting effects would be sup- 
pressed (there are no observable implications for splitting of 
isotropically emitting sources). 

4. CONCLUSIONS 

In this paper we have described the effect of photon-ALP 
beam splitting in the presence of an external magnetic field - 
an effect which arises solely by virtue of the interaction of the 
particle with the EM field. This effect is very different than 
(but physically related to) that of photon-particle oscillations 
which were discussed in C09, and provides a complementary 
means for probing axion and ALP physics in magnetized envi- 
ronments. In this paper we focus on the observable signatures 
of this effect in the light curve of compact, highly magnetized 
objects and, in particular, the case of magnetars. The formal- 
ism developed here, and applied to magnetars, is limited to 
the case of negligible photon and ALP mass (relative to the 
interaction term) or at resonance where they equate. 



We have shown that, by studying the radio light-curves of 
magnetars, one can be sensitive to light bosons down to very 
low values of the coupling constant - about 3-4 orders of mag- 
nitude smaller than current CAST limits and other indirect as- 
trophysical constraints. The parameter space probed by such 
experiments is shown in figure[3][for various combinations of 
£ = (foB /10 16 G)(A/m)] and nicely complements spectro- 
scopic searches for photon-particle oscillation features (C09). 
Note that, for a given £, there is a maximum g beyond which 
one beam is attenuated (equation QT| is not satisfied). In this 
case, one can probe pulse shifting, as discussed above. 

The detection and verification of the splitting and shifting 
effects may be done by looking for a typical phase differences 
between pulses in the light curves of magnetars. In particu- 
lar, the splitting/shifting phase would be wavelength depen- 
dent and would diminish at shorter wavelengths. In addition, 
the fluxes and polarizations of the split pulses would be sim- 
ilar. Split or shifted pulses are likely to have a different po- 
larization than pulses whose photons do not mix with bosons. 
Splitting and shifting effects could also depend on the rota- 
tion phase of the magnetar reflecting, perhaps, the configu- 
ration of the magnetic field along the line-of-sight. The var- 
ious trends discussed above provide simple yet robust tests 
for photon-particle mixing. It should be noted, however, that 
failing to detect beam splitting and shifting cannot be used to 
confidently exclude part of the axion/ALP parameter space so 
long as our understanding of the radio emission from magne- 
tars is incomplete. Finally, we wish to note that the results 
presented here are general and apply to all (astrophysical) set- 
tings in which beamed, long wavelength emission originates 
from highly magnetized regions. 
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5. APPENDIX 

Consider the scalar QED Lagrangian (neglecting mass terms and magnetic fields; Guendelman 2008a,b): 

L = y/*d » y/ + ieA ( y/* d y/ ~ y/<9 y/* ) (12) 

e is the charge of the scalar particle and A the electric potential. Here \\r = [\a) + i\y)]/V2. Using Euler-Lagrange equations we 
find, 

□ y/-2/M°<? y/ = 0^ (k 2 - co 2 - 2ea>A°) y/ = (13) 

where the last step used an ansatz such that y/ = e~'°"yf(x). This expression is completely analogous to that given in equation 
[3](after diagonalization) with eA° — > gB\\/2. The eigenstates defined in equation [5] correspond to different charge states in the 
scalar QED picture and an analogue to the electric field exerts an opposite force on each charge state leading to the beam splitting 
effect shown in figure 1 . 



